Landau-Zener problem for energies close to potential crossing points. 
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■ We examine one overlooked in previous investigations aspect of well - known Landau - Zener (LZ) 

O ' problem, namely, the behavior in the intermediate, i.e. close to a crossing point, energy region, 

when all four LZ states are coupled and should be taken into account. We calculate the 4x4 
connection matrix in this intermediate energy region, possessing the same block structure as the 
known connection matrices for the tunneling and in the over-barrier regions of the energy, and 
"j^ ' continously matching those in the corresponding energy regions. Applications of the results may 

concern the various systems of physics, chemistry or biology, ranging from molecular magnets and 
glasses to Bose condensed atomic gases. 
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' ^ ' Standard textbook LZ theory Q treats of two linear diabatic potentials U'^ ± FX crossing problem {X = is the 
Q : crossing point). However, in spite of more than half of century history, semiclassical solutions of this problem have 
^ ' been found only in the limits of small or large energies E (wc will term these regions as tunneling or over-barrier, 
respectively), i.e. for 



|C/#-£;|»C/i2; |C/#-S|«C/i2, (1) 
where U12 is inter level interaction which in LZ model does not depend on X. For the intermediate energy region 



> 

^ ■ \U*-E\< U12 (2) 



there known only interpolating relations between exponentially decaying solutions in the tunneling and oscillating 
solutions in the over-barrier energy regions (see e.g. 0). Analytical and numerical study of this region Q is the 



^ . objective of this paper. 

Our approach is motivated by semiclassical instanton approximation , Q , 0] . The idea is to construct two linearly 
^ , independent continuous (with continuous first derivatives) approximate solutions to the Schrodinger equation, which 
in the asymptotic region coincide with semiclassic solutions, and in the vicinity of the turning points - with the exact 
' ^ , solutions of the so-called comparison equation (i.e. the exact solution of the Schrodinger equation for the chosen 
appropriately approximate near the turning points potentials). In what follows the Weber equation jOj will be used 
^ , as the comparison equation, valid in the vicinity of the second order turning points for an anharmonic potential 0] , 
0,0. To justify this choice it is sufficient to note that anharmonic corrections remain semiclassically small (i.e., 
proportional to higher orders of h series) in the region where the solutions of the comparison equation have to be 
matched smoothly with the semiclassical solutions. Luckily the analogous approach is valid to treat two diabatic 
potential crossing point (LZ problem), and the comparison equations for this case are two coupled Weber equations 
with the indices and arguments determined by the solutions of algebraic characteristic equation. 

LZ problem for crossing diabatic potentials is equivalent to the coupled Schrodinger equations which can be trans- 
formed by the substitution 

* = exp(KX)4' (3) 
into the 4-th order linear differential equation with independent of coordinates coefficients at the derivatives 

+ 4k£>^$ + {6k^ - 207^)1)2$ + 4:{k^ - a-f'^n - ^j^f)D<^+ (4) 
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2a7 K - 27^/K + 7'*(a^ 



'■12 



where D" = and 7 ^ 1 is the dimensionless semiclassical parameter which is determined by the ratio of 

the characteristic potential scale over the zero oscillation energy, and all other dimensionless appropriately rescaled 
variables are 



U*~E ^ 2aoF 



2U^ 
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(5) 



where scale for the energy is given by fl^ — F'^/mUu {m is a mass), and space scale is determined by the characteristic 
size oq of the potential in the vicinity of the crossing point. 

The equation admits semiclassical solutions by Fedoryuk method ^3 " ^3 since the coefficients at the n-th 
order derivatives proportional to 7^", and all the four asymptotic solutions read as 



(sc) 



Xj{x)dx ] ,j = (++,+-,-+, ), 



(6) 



where we designated Xj = A" + Uj, A° = ± ^'^^(0^"^/^!^^^^^^^, and uj = 7/((A°)^ — aj) ^. 

The equation Q, up to anharmonic terms proportional to X^D^, can be formally derived by simple 

manipulations (two sequential differentiations and summations) from the following second order equation 



where the coefficients are 



ao = K 



D"^^ + (ao + aiX + aaX^)^ = , 



2k 



where the equation for k referred in what follows by the characteristic equation is 

- aj^K^ + \l^u\^ = -k4(52(1 + 2(5) + R{k, 8) , 

where 



(7) 



(8) 



(9) 



and 



i?(K, b) = (2kS)-1(1 - 3J)(1 + -Q- v/l-2Q2)) ; Q 



The fundamental solutions to Q read as 



± 



1/4 



X 



1 

2k 



where 



^4^2 X -1/2 



ao 



4a2 



(10) 



(11) 



(12) 



In the tunneling and over-barrier regions of energies, where 5 < 1/4, these 4 solutions (2 solutions of pi|l for two 
largest modulus roots of the characteristic equation ^) can be separated into two independent pairs (orthogonality 
of the Weber functions with different indices). Thus one can say that the crossing point is equivalent to two second 
order turning points with different Stokes constants (see e.g., |l3l|1. 

In the tunneling region the two largest modulus roots of are (two other roots are small and do not satisfy 
semiclassical approach) 



K — ±Ko 1 ± 



2 2k^ - 072 



1/2 



(13) 
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and the four linearly independent solutions 1)11(1 are matched to the semiclassical solutions ((Sj) in the region a > 
f\X\ > Mi2, where the exponent of the wave function can be expanded over small parameter S 

$ cx exp (^kX + S{KXf - ^S^inXf + ^ . (14) 

Since k\X\ < {4S)^^, the convergent with alternating signs expansion (|13|l determines the accuracy of the asymp- 
totically smooth transformation. Putting altogether we end up at the conclusion that anharmonic corrections to the 
Weber functions ((ll|l are small (by other words the parameter 5 determines the accuracy of our approximation) . The 
same kind of analysis can be performed in the over-barrier region, where one finds two imaginary largest modulus 
roots of the characteristic equation (see details in ^S] ) . 

More difficult task is to find solutions in the intermediate energy region, where two roots of the characteristic 
equation are real and two are imaginary ones having the same modulus, i.e. moving upon a variation along a circle 
with the radius 7-^/ U12/2. In this case the semiclassical solutions can be presented as certain linear combinations of 
the comparison equation solutions, and the roots are 

Ki,2 ^ ±ld — exp(iip) ; K3,4 ^ ±i"fd — e-Kp[-iip) , (15) 

where 



tamp — \ . (16) 

^ \ ui2 + a ^ ^ 

Correspondingly to these roots (|15|l the arguments and the indices of the Weber functions (|ll|l . (|12|) read as 

zi = Z2 = 2Kint\/%^texp{-i(p/2){X - {2Kinty^ exp{-iip) ; (17) 



Z3 = -Z4 = 2K™t y^Sint exp{i(p/2){X - {2Kint) ^ exp(z(/9) , 



and 



Pi = P2 - 1 = -1 



■ exp(-i(p)(l 4- 2(5f„t exp(-2i(/?)) ; 



(18) 



P4 = P3 - 1 



■ exp(i<^)(l -f 25'f^^ exp(2i(^)) , 



where 



(19) 



The semiclassical solutions © are matched asymptotically smoothly to the linear combinations of the Weber functions 
in the region ui2 > f\X\. Since from l)17|l . 118|) follow that at large Sint the indices of the Weber functions are also 
large, one can use known due to Olver 14], 15] asymptotics of the Weber functions with large arguments and indices 



Dp{z) cx exp 



A p 



1/2 



dz 



(20) 



At z'^ 3> 4]p + (1/2)] H20() is reduced to the usual asymptotic expansion of the large argument Weber functions, and 
in the opposite limit (i.e. in the intermediate region) H2()|l corresponds to the expansion of the exponent over odd 
powers of z. We can also find asymptotics to the solutions of lO 



$0 oc exp ^— i J ^/ao + aiX + a2X'^dx 



(21) 



valid at arbitrary values of the parameters Ui (02 — including). 
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The fundamental solutions to the comparison equation are the asymptotics for the wave functions in the form 
namely, 



(22) 



and 



cx exp(Fi(X)) , *3 cc eM~FiiX)) , ^T, *3 « eM^F^iX)) , cx exp(-zF2(X)) , 



where 



Fi,2(^) = 7V"i2±a(l + (5„t)X - Kf^tSf^t exp{~2iip)X' + 



12U12VW12 ± a V ""12 



(23) 



(24) 



The wave functions (|23l) asymptotically smoothly turn into the semiclassical functions (0. The accuracy of this 
matching is determined by the anharmonic corrections, i.e. by the parameter dint (|19l) . and the Olver asymptotic l|2U|) 
works even on the boundary ~ 4\p + (1/2) |. The parameter Smt is no more small one, when simultaneously 



lal < 



2/3 



Ul2 
272 



1/3 



U12 < 



(25) 



However the asymptotic matching of the solutions should be performed at small \X\ < 7"^, where the comparison 
equation p2(l . and, therefore, the characteristic equation ^ are valid, even though upon increasing Smt the potential 
becomes more and more anharmonic one. At a = 0, and M12 — 0, the equation ^ (with the term R{k, S) taking into 
account) has the doubly degenerate root k = 0, i.e. in terms of H12|l ai = ±7'^/, 02 = 0. Thus in this limit (|12|l is 
equivalent to two decoupled Airy equations, corresponding to the diabatic potentials. These solutions turns smoothly 
into the semiclassical ones (O, and the anharmonic corrections in the matching region are small over the parameter 
7-1/2. 

We conclude that in the both intermediate energy subregions: large k (i.e. cx 7), and small small k (i.e. cx ^/7) 
the comparison equation (|12(l is reduced to two decoupled equations, Weber or Airy ones, respectively. This simple 
observation enables us to construct the universal connection matrix for the both intermediate energy region by using 
Olver asymptotic expansion (|21|) . The four roots H15|) distributed over the circle with the radius 7-^/ (M12/2) on the 
comples plane determine the following combinations of the comparison equation solutions matching the semiclassical 
solutions jnj. Namely 



-+ ' 



*7 + *3 



(26) 



where the superscript sc points out the semiclassical solutions. Combining together the asymptotic expansions for 
these combinations, we find at the crossing point, the matrix [/" is 



U'' 



(V2Vr(g*))cxp(-2x(<z*)) 



exp(— 27r(72) cos(7rgi) 



(r(g)/^ 





/27r) exp(2x((7))(l - exp(-27rq2) coa^{TTqi)) 
exp(— 27rq2) cos(7rgi) 




(27) 



where 



— exp(— 27rg2) cos(7rgi) 

— exp(— 27r(72) cos(7rgi) 

(V2^/r(g))exp(2x(<z)) 

{T{q*)/V2^) exp(2x(g*))(l - cxp(~27rg2) cos2(7rgi)) J 



7U12VW12 ± a , 
q = qi+iq2] gi,2 = ; g = gi - «g2 , 



(28) 



and, besides, we introduce the following abridged notations 



X = xi + «X2 ; 2x1 = gi - ( gi - 2 ) l^l + 'PI2 ■. 



(29) 
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and analogously 



2X2 = (72 - 92 In |g| - V? ( gi - i ) , (30) 



where (p is defined by 

The connection matrix (|27() in the intermediate energy rerion is our main result and the motivation of this pub- 
lication. This matrix generalizes the results we presented in It is ready for further applications and to reap the 
fruits of the result we compute the LZ transition probability |Tp universally valid for the tunneling, over-barrier and 
intermediate energy regions (solid line in the Fig. 1). It is instructive to compare our result with the perturbative 
Landau approach (see e.g., valid at small coupling constants. In the first order perturbation theory the transition 
amplitude A^^^ reads 

4'l = 2^^il2y^^^(22/3e), (31) 

where we designated e — — a/^^/"^, and Ai is the first kind Airy function. All even higher order terms equal zero, and 
odd terms read as 

A^ir'^ = (4'i)'"+^ ■ (32) 
Therefore the series can be easily sum up giving the generalized Landau formula 

ALz = A'^ll[l + {ALzfY^. (33) 

We show this perturbative solution by the dashed line on the Fig. 1. Note that although the equation (|32(l reproduces 
the oscillating energy dependence of the LZ transition amplitudes, the equation gives the period of the oscillations 
which are quite different from those we calculated by our 4x4 connection matrix H27|) . The difference occurs because 
the perturbation (say 2) method disregards the contributions from the increasing solutions to the Schrodinger equation, 
which are relevant in the intermediate energy region. 

The research described in this publication was made possible in part by RFFR Grants. One of us (E.K.) is 
indebted to INTAS Grant (under No. 01-0105) for partial support, and V.B. and E.V. are thankful to CRDF Grant 
RU-C1-2575-MO-04. 
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Figure Captions. 

Fig. 1 

Energy dependent LZ transition probability: 

Solid line - the 4x4 connection matrix H27|) calculations; 

Dashed line - the generalized perturbative Landau formula H33|l . 




Fig.l 



